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ABSTRACT 



One of the generic predictions of modern cosmological 
models is that large galaxies should have experienced 
many mergers with smaller galaxies at some point in their 
past. Debris from such encounters will leave spatially 
distinct substructure in the stellar haloes of nearby galax- 
ies, detectable for a few orbital periods after the final 
merger. In the case of the Milky Way, kinematic data 
from surveys such as RAVE and satellite missions such 
as GAIA will allow us to probe much more of the merger 
history, and to connect the properties of the stellar halo 
with those of local dwarf galaxies. To estimate what 
these programmes may discover, we review current ob- 
servations of minor mergers in nearby galaxies, and com- 
pare these with predictions from a semi-analytic model 
of galaxy formation. 
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1. INTRODUCTION: THE ROLE OF MERGERS 
IN CURRENT COSMOLOGICAL MODELS 



In cosmological models consistent with current observa- 
tions of the cosmic microwave background (CMB), large- 
scale structure, and estimates of the mean baryon density, 
roughly 85% the matter in the universe is non-baryonic 
cold dark matter (CDM). Small fluctuations in the ini- 
tial distribution of CDM are known to lead, via gravita- 
tional instability, to the formation of dense, gravitation- 
ally bound structures at later times. These dark matter 
'haloes' should reach masses and densities sufficient to 
allow gas to cool and condense within them, thereby per- 
mitting efficient and sustained star formation, at redshifts 
or ^20 or higher. Thus they represent a natural structural 
framework within which visible galaxies can form, from 
the earliest epoch onwards. 

A generic prediction of cosmological models based on 
CDM is that the dark matter haloes around galaxies con- 
tinue to grow through accretion and mergers with smaller 
haloes, right up to the present day. The merger rate 
for galaxy haloes reaches a peak at a redshift of z ~1- 
3, however, and more recent mergers between galaxy 



haloes and their associated galaxies will typically be mi- 
nor ones, adding little to the total mass of the larger com- 
ponent. The average merger in the local universe may 
consist of a small, low surface-brightness dwarf galaxy 
being stretched out and disrupted over the course of sev- 
eral orbits around its parent. The only direct signature of 
such an event would be minor fluctuations in the surface 
brightness of the stellar halo around the larger galaxy, 
lasting for a few orbital periods (Johnston et al. 2001). 
Thus, while low-redshift mergers provide an important 
empirical test of the CDM framework, they may be ex- 
tremely difficult to detect in practice. 

In the case of the Milky Way, stellar kinematics pro- 
vide a unique opportunity to study the merger history in 
much greater detail, identifying debris from older or more 
minor mergers thanks to its clustering in phase-space 
(Tremaine 1993; Johnston 1998; Helmi & White 1999). 
Ground-based radial velocity surveys such as RAVE 1 and 
future astrometric missions such as GAIA 2 will explore 
a completely new range of parameter space in the merger 
history of this one particular galaxy, and should thereby 
provide an important link between local stellar popula- 
tions and the properties of small galaxies forming at the 
highest redshifts. To explore the potential of these pro- 
grammes, we will review current observations of debris 
from minor mergers in external galaxies, and attempt to 
extend these results to older or smaller satellites using a 
semi-analytic model of galaxy formation. 



2. OBSERVED STREAMS IN NEARBY GALAX- 
IES 



Given their theoretical significance, there has been much 
recent interest in detecting ongoing minor mergers in lo- 
cal galaxies. The most spectacular discoveries of this 
kind were those made in our own galaxy and in M31. 
In 1994, Ibata et al. announced the discovery of the 
Sagittarius Dwarf, a small system in the process of be- 
ing disrupted by the Milky Way, whose tidal debris is 
now spread out over a large fraction of the sky. More 
recently, and equally spectacular tidal stream was dis- 
covered around M31 (Ibata et al. 2001; McConnachie 

1 http://www.aip.de/RAVE 
2 http://as tro.estec.esa.nl/GAIA 
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et al. 2003), and there is growing evidence for another 
stream in the plane of the Milky Way (e.g. Newberg et 
al. 2002; Ibata et al. 2003). The LMC-SMC system can 
also be considered an ongoing minor merger, since it has 
already experienced strong tidal effects due to its prox- 
imity to the Milky Way (e.g. Harris & Zaritsky 2004), 
and will eventually merge with the Galaxy as dynamical 
friction acts on its orbit. 

Much of the tidal debris detected locally has a very low 
surface brightness, and would be difficult or impossi- 
ble to detect in distant systems. Nonetheless, there has 
been some progress in detecting debris around galax- 
ies outside the Local Group. Early work by Malin & 
Hadly (1997) found faint features in the stellar haloes of 
several galaxies, including M83, M104 and NGC2855, 
some or all of which may be tidal debris from minor 
mergers. A clear example of a tidal stream was discov- 
ered around the edge-on spiral NGC5907 by Shang et 
al. (1998), and the nearby active galaxy Centaurus A has 
a tidal feature of young blue stars, possibly associated 
with the recent merger that triggered its central activity 
(Peng et al. 2002). These results were reviewed and dis- 
cussed recently by Pohlen et al. (2004), who also report 
newly detected streams in four other galaxies. Their sys- 
tematic search around 80 edge-on disk galaxies yielded 
only one of these examples, however, suggesting the in- 
cidence of obvious streams is low. More recently, Forbes 
et al. (2004) have reported the serendipitous discovery of 
a dwarf in the process of disruption, with well-defined 
tidal streams, in the background of a Hubble Advanced 
Camera image of the Tadpole system. These discover- 
ies confirm that minor mergers are still taking place, but 
suggest that it may be hard to obtain a large sample to 
compare with theoretical models. Finally, we note that 
stellar streams are not the only tracer of minor mergers in 
external galaxies; others include distinct populations of 
globular clusters, gas or dust lanes in early-type galaxies, 
and kinematic or structural anomalies. The theoretical in- 
terpretation of these features is less straightforward, how- 
ever, as their properties may depend more on the long- 
term response of the main system (e.g. via AGN or star- 
bursts) than on the nature of the initial perturber. 



3. THE UNRESOLVED PROBLEM OF DWARF 
GALAXY FORMATION 



Given this growing body of information on minor merg- 
ers and tidal streams, it is worth estimating how often 
these events are expected to occur in current galaxy for- 
mation models. In recent years several groups have begun 
to study the formation of the stellar halo in its full cos- 
mological context (e.g. Bullock, Kravtsov, & Weinberg 
2001; Bekki & Chiba 2001; Brook et al. 2004; Bullock 
& Johnston 2004). The problem is not as straightforward 
as it seems, however. The properties of the stellar halo 
will depend on when and how the satellites that merged 
into it first formed. While simulations of structure forma- 
tion have now converged on the properties of dark matter 
structures on the relevant scales, it is still not clear how 
these structures are populated with visible stars. 

It has been known since the early days of CDM that 



galaxy formation must be systematically less efficient 
in smaller haloes, since the luminosity function of field 
galaxies flattens below some magnitude, whereas the 
halo mass function should be close to a power law over 
the corresponding range (e.g. Kauffmann et al. 1993). 
More recently, high-resolution simulations by Klypin et 
al. (1999) and Moore et al. (1999) have demonstrated 
that this discrepancy is even more dramatic for satellites 
within systems like the Local Group. Using internal ve- 
locities to relate dark matter subhaloes to visible systems, 
they showed that down to the scale of the smallest sys- 
tems, there should be almost 100 times more subhaloes 
around the Milky Way than there are detected satellites. 
Subsequently, there has been much discussion of whether 
dwarf galaxies populate a small subset of all subhaloes 
over a wide range of mass/internal circular velocity, as 
proposed by Moore et al. (1999), or whether they corre- 
spond to kinematically cold stellar cores within the most 
massive subhaloes, as suggested by Stoehr et al. (2002) 
and Hayashi et al. (2003). 

The second of these solutions is appealing, as it could 
indicate that there is a single characteristic halo mass be- 
low which star formation ceases. This solution is strongly 
disfavoured, however, if not ruled out, by the strong spa- 
tial clustering of the satellites of the Milky Way. The 
best estimates of the mass of the Milky Way put its virial 
radius somewhere around 300 kpc. If this is the case, 
then two-thirds of its satellites lie within the central third 
of the halo. The chance of this happening if satellites 
populated only the most massive haloes is less than 1% 
(Taylor et al. 2003). The disagreement in the cumulative 
radial distributions is shown in Figure ^ where we com- 
pare the average results for the dozen most massive satel- 
lites from each of a large set of semi-analytic models of 
halo substructure, generated using the method outlined in 
Taylor & Babul (2004a, 2004b, 2004c) (solid line + dot- 
ted 99% contours) with the distribution for the satellites 
of the Milky Way (upper line & triangles). Thus at least 
some of the satellites must reside in smaller haloes. 
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Figure 1. The cumulative radial distribution of visible 
satellites around the Milky Way (upper line/triangles), 
compared with the predicted distribution of the most mas- 
sive subhaloes (lower solid line, with dotted 99% con- 
tours. 

There may be additional clues to the origin of the local 
dwarfs in their kinematics, which appear very different 
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from those of randomly selected subhaloes. In Figure[2] 
for instance, we show how the positions and radial veloc- 
ities of known satellites (points with error bars) compared 
with the distribution of massive haloes taken from a large 
set of semi-analytic haloes. There is a clear mis-match 
between the two distributions in projection along either 
axis. 
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Figure 2. Radial velocity versus position for the satellites 
of the Milky Way (points with error bars) and the most 
massive subhaloes from a series of semi- analytic haloes 
( smaller symbols). 



One possible explanation, both for the spatial clustering 
and for the observed velocity distribution, is that most 
of the dwarfs are old, and that the lower cutoff to dwarf 
formation increases with time. In particular, Kravtsov et 
al. (2004) have produced a detailed model of dwarf for- 
mation where this occurs naturally, through a combina- 
tion of internal and external feedback. In the next section 
we will consider a simplified version of this model and 
discuss its implications for stream formation. 



A SEMI-ANALYTIC MODEL OF STREAM 
FORMATION 



number and spatial distribution of surviving satellites at 
the present-day. 

Having identified which systems form stars, we can then 
estimate when stellar material will be stripped from them 
and added to the stellar halo. We assume the stars within 
a given subhalo are restricted to its central region, in 
keeping with the observed sizes of dwarf galaxies, and 
consequently that tidal stripping does not produce visi- 
ble streams until a system has lost most of its original 
dark matter. Considering the stellar-to-total mass ratios 
of larger galaxies, we set a mass-loss threshold at 83%, 
and refer to systems that have lost less than this fraction 
of their original mass as 'surviving satellites', and those 
that have lost more as 'tidal streams'. Systems that lose 
more than 99% of their mass may become completely un- 
bound (Taylor & Babul 2004a), so we refer to them as 
'disrupted' in what follows. 

The distribution of stripped material from systems con- 
taining stars is shown in Figure[5] for one particular semi- 
analytic halo. The debris is plotted at the point where 
it was first stripped off, to highlight the underlying or- 
bital structure. Over time this material will be heated 
and phase-mixed, producing a smoother final halo. The 
colours in the figure indicate the age at which the stripped 
system first merged with the main halo. There is a strong 
age gradient in the final system, with ancient mergers 
contributing most of the material within the solar radius, 
while more recent mergers add streams to the outer halo. 




Given the complexity of modelling dwarf galaxy forma- 
tion self-consistently, it is useful to construct a simpli- 
fied semi-analytic model of the process, to get a prelim- 
inary estimate of how often minor mergers with dwarf 
satellites would have occurred in a system like the Milky 
Way. Our model is based on the semi-analytic model of 
halo formation presented in (Taylor & Babul 2004a,b,c), 
which predicts the numbers, orbits, internal properties 
and dynamical evolution of dark matter subhaloes within 
a galaxy halo. To make predictions about visible satel- 
lites, we will suppose that stars form preferentially in sub- 
haloes with deeper potential wells, but also that the lower 
threshold to this process increases with time, such that 
the some of the largest subhaloes around the Milky Way 
have never formed stars. Specifically, we choose a time- 
varying threshold in peak circular velocity of the form 
V p > Vp.o (1 + z)~ a . Subhaloes that exceed this thresh- 
old at any time before they merge with the main halo are 
assumed to host visible dwarf galaxies. The two param- 
eters in this model can be adjusted to give the right total 



Figure 3. Tidal debris from merging satellites in a semi- 
analytic model halo, plotted at the point where it was first 
stripped. The material is colour-coded by age, the central 
material being systematically older. The figure is in the 
plane of the disk and the scale is in kpc. 



5. RESULTS: PREDICTED STREAM PROPER- 
TIES 



Using the model outlined above, we have generated 
streams in a large set of model haloes. We find that we 
can produce an average of ten 'surviving' satellites per 
halo, roughly the number known to-date for the Milky 
Way, by setting U Pi o = 135 kms -1 and a — 2/3 in 



4 



the expression given above. With this choice of param- 
eters, our model predicts ^8 'tidal streams' (i.e. heav- 
ily stripped stellar systems) per halo, and more than 100 
disrupted dwarf galaxies, as well as 15 that have merged 
directly into the centre of the main galaxy on very ra- 
dial orbits. Figure |4] shows the distributions of various 
parameters, including the original total mass, formation 
redshift, merger redshift, and final radius (four columns, 
from left to right) for these four classes of objects (four 
rows, from top to bottom). Various interesting trends ap- 
pear in these distributions; in particular, the progenitors 
of tidal streams (second row) lie at somewhat smaller 
radii (~ 30 kpc), are slightly older, and have larger origi- 
nal masses than the surviving satellites (top row). 
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Figure 4. Distribution of original total mass, forma- 
tion redshift, merger redshift, and final radial position 
( columns from left to right) of stellar systems that survive 
at the present-day, tidal stream progenitors, disrupted 
systems and central mergers (four rows, from top to bot- 
tom) in a large set of semi-analytic haloes. 

We can compare these predictions with the observations 
discussed above. On the one hand, only a few streams 
have been detected in the Milky Way (Sagittarius, Canis, 
and possibly a few older streams associated with objects 
such as the massive globular cluster uj Cen; cf. Bekki & 
Freeman 2003). On the other hand, given the difficulty 
of detecting older streams, substantial incompleteness is 
not implausible. Our results suggest that GAIA may ex- 
pect to sample stars from a much larger population of 
disrupted systems - on the order of 100. Most of these 
systems will have merged with the halo at redshifts of 2 
or higher, however, and will consist of old stars, strongly 
mixed in phase-space by repeated scattering from their 
initial orbits. It will take much more detailed modelling 
of the disruption and mixing process (cf. Johnston 1998; 
Helmi & White 1999) to determine how many distinct 
structures should be detectable in practice. 
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